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ABSTRACT 



By using bidimensional spectral data obtained at the 6m telescope for the 
Virgo spirals NGC 4216 and NGC 4501, we have found chemically distinct 
metal-rich nuclei in these galaxies. Under the assumption of equal ages for 
the nuclear and bulge stellar populations, the metallicity difference between 
the nuclei and their environments in the galaxies is estimated as a factor of 2. 
But we have also found an age difference between the nucleus and the bulge 
in NGC 4216: age-metaUicity disentangling on the diagrams (H/3, Mgb), (H/9, 
[MgFe]), and (H/?, < Fe >) results in an age estimate for the nucleus of 8-12 
billion years, the bulge being older by a factor of 1.5-2; and the self-consistent 
metallicity difference estimate is then a factor of 3. The solar magnesium-to-iron 
ratios in the galactic nuclei show evidence for a long duration of the secondary 
nuclear star formation bursts which produced the chemically distinct stellar 
subsystems. Detailed morphological and kinematical analyses made for the 
stellar and gaseous structures in the centers of NGC 4216 and 4501 have 
revealed the presence of circummiclcar stellar-gaseous disks with radius of some 
hundreds parsecs which demonstrate fast axisymmetric rotation and lie exactly 
in the planes of the main galactic disks. 



Subject headings: galaxies: spiral — galaxies: individual (NGC 4216) — 
galaxies: individual (NGC 4501) — galaxies: abundances — galaxies: nuclei — 
galaxies: kinematics and dynamics — galaxies: structure 
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1. Introduction 

We search for chemically distinct stellar nuclei in disk galaxies for several years, and 
this search has been successful: chemically distinct nuclei are found in a dozen of disk 
galaxies (ISiPchenko 19971 ) . From the beginning we tried to find a relation between the 



presence of chemically distinct nuclei and the galactic structure. At first we have noticed 
that isophote twist is absent, as a rule, in the galaxies with chemically distinct nuclei 
( [Sirchenko fc Vlasyuk 1992| ); so we have associated the chemically distinct nuclei with 
central mass concentrations within axisymmetric bulges and have thought them to be the 
result of the internal evolution. But later some arguments have been accumulated that the 
presence of chemically decoupled nucleus is often accompanied by kinematically decoupled 
gaseous subsystems, such as polar rings, so it may be a consequence of a past interaction 
with some gas accretion and subsequent secondary star formation burst in the nucleus 
( ISirchenko et al. 1997| , Fil'chenko 19981) . 



The chemically distinct nuclei in elliptical galaxies are more widely known, and their 
presence is often related to nuclear compact stellar disks within main galactic bodies. 
Bender and Surma ( |1992| ) reported the discovery of a magnesium index enhancement in 
dynamically decoupled cores of four elliptical galaxies. Since they deal with dynamically 
hot elliptical galaxies, the chemically decoupled cores have been easily identified as being 
disks from dynamical arguments, namely, from their fast rotation and, in two cases, from 
a slight decrease of stellar velocity dispersion toward the centers. Later for two galaxies a 
photometric study has confirmed the presence of the nuclear stellar disks with appropriate 
sizes (NGC 5322: [Scorza fc Bender 1995| ; NGC 4365: [Surma fc Bender 1995| ). Recently 



another two elliptical galaxies, NGC 4816 and IC 4051 in the Coma cluster, were found to 
possess chemically (and kinematically) distinct cores which are in fact central stellar disks 
( [Mehlert et al. 1998|) . The hypoth esis that the kinematically and chemically distinct cores 
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in elliptical galaxies are nuclear stellar disks formed via dissipational major merger events 
has become popular. However one must admit that there exists alternative points of view: 
Carollo et al. (|1997|) who used the HST WPFC2 data have not found any nuclear disks 
in the elliptical galaxies with kinematically and chemically distinct cores NGC 2434 and 
NGC 7192 and have concluded that magnesium-rich cores "may arise from a large variety 
of morphological and dynamical structures" . 

Spiral galaxies were studied less thoroughly, but we know one example which was 
considered in detail. It is NGC 4594, which has been investigated both spectroscopically 
and photometrically more than once. Particularly, results of bidimensional spectroscopy 
presented by Emsellem et al. ( 1996 ) are of great interest for us in the frame of our topic. 



NGC 4594 has an inner disk with the radius of 11" which is disconnected to the more outer 
main disk of the galaxy - this was already shown by Burkhead ( |1991| ). Two-dimensional 
maps of absorption-line indices Mgb, Fe5270, Fe5335, and Fe5406 obtained by Emsellem 
et al. (|1996|) have revealed the presence of a chemically distinct nucleus in NGC 4594; 
but, interestingly, the magnesium- index difference between the nucleus and the inner disk 
is smaller by a factor of 2 than the corresponding difference between the nucleus and the 
bulge (along the minor axis), and there are no difference of iron indices between the nucleus 
and the inner disk. In some sense, just the inner disk in NGC 4594 represents a chemically 
distinct core. 

So, we suspect that the presence of chemically distinct nuclei in spiral galaxies may 
be related to their inner stellar disks; may be, it is the same thing. But the sample of 
spiral galaxies with decoupled circumnuclear stellar disks is not larger (perhaps, even 
smaller) than the sample of spiral galaxies with chemically distinct nuclei. An alternative 



starting point is proposed by a recent work of Rubin et al. ( |1997| ) who have reported 



the discovery of 14 rapidly rotating, discrete circumnuclear gaseous disks in Virgo cluster 
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spiral galaxies. Among those, NGC 4216 possesses one of the most extended solid-body 
rotating circumnuclear gaseous disks: Rubin et al. (|1997|) give its radius equal to 10" (or 
800 pc) and its angular rotation velocity equal to 28 km/s/arcsec, or 350 km/s/kpc. Since 
NGC 4216 is also included into our list of candidates for possessing a chemically distinct 
nucleus (ISiPchenko 1994|) , we think it to be an appropriate object which may contain a 
distinct nuclear stellar disk; due to the large extension of the circumnuclear gaseous disk, 
the stellar disk may be also large enough and so may be resolved under moderate seeing 
conditions. Another Virgo spiral, NGC 4501, is not mentioned by Rubin et al. ( 1997 ) as 
having a circumnuclear gaseous disk; moreover. Keel ( [L983|) stated that its central emission 
is confined to the nucleus and does not show any disk structure. But it is also in our 
list of candidates for possessing a chemically distinct nucleus ( ISiPchenko 1994| ), and we 
have decided to study it more carefully together with NGC 4216. We have undertaken 
bidimensional spectroscopy of the central parts of these galaxies. Section 2 describes 
the details of our spectral observations and of other observational data involved into our 
analysis. In Section 3 we discuss the morphology of surface brightness distributions in 
NGC 4216 and NGC 4501. Section 4 deals with the kinematics of the gas and stars in the 
central parts of these galaxies. Section 5 presents the discovery of their chemically distinct 
nuclei. Finally, Section 6 contains a brief summary and conclusions. 



2. Observations and Data Reductions 

We have undertaken spectral investigation of NGC 4216 and NGC 4501 at the 6m 
telescope of the Special Astrophysical Observatory of the Russian Academy of Sciences. A 
full journal of the spectral observations is presented in Table 1. 

In July 1996 we observed NGC 4501 with the Long-Slit Spectrograph of the 6m 
telescope, under rather good seeing conditions. The slit, with a length of 2'5 and a width 
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of 2", has been set in three position angles: along the major axis, along the minor axis, and 
between them, at 45° from the major axis. We exposed a red spectral range containing 
strong emission lines Ha and [N11]A6583 and derived line-of-sight velocities of the ionized 
gas along the slit by measuring barycenters of these emission lines. The dispersion was 
1.56 A per pixel, and the spectral resolution varied along the slit between 3.2 and 3.9 A. 
Test measurements of the sky emission line [OI]A6300 have shown that the accuracy of our 
line-of-sight velocity measurements is 6-8 km/s. In the center where the emission lines are 
especially strong we have used one-pixel spatial bin equal to 0'.'4, in more distant regions - 
three-pixel binning. 

The more detailed spectral investigation of the innermost parts of NGC 4216 and 
NGC 4501 has been undertaken with the Multi-Pupil Field Spectrograph (MPFS) of the 
6m telescope (for description of the early variant of the MPFS - see [Afanasiev et al. 19901 ). 
This type of spectrographs is still rather untraditional; a detailed explanation of all the 
profits of bidimensional spectroscopy can be found in the paper of Bacon et al. (|1995|) , 



and some details of our particular approach are given in our work ( ^il'chenko et al. 1997| ). 
Briefly, we obtain simultaneously about one hundred spectra, each spectrum being from a 
square spatial element of 1'.'3 x 1'.'3. The total element array covers a rectangular area in the 
center of a galaxy. The basic data reduction steps - bias subtraction, fiatfielding, cosmic 
ray hit removing, extraction of one-dimensional spectra, transformation into wavelength 
scale, construction of surface brightness maps and velocity fields - are performed using the 
software developed in the Special Astrophysical Observatory ( [Vlasyuk 1993|) . Observing 
NGC 4216 and NGC 4501, we used two different setups. In May 1996 and in January 1998 
a CCD of 520 x 580 was in work; we exposed an area of 8 x 12 elements (10" x 16") with a 
dispersion of 1.6-1.7 A per pixel (spectral resolution being 3-5 A). In May 1997 a CCD of 
1040 X 1160 was operated, and we exposed an area of 8 x 16 elements (10" x 21") with a 
dispersion of 0.9 A per pixel and a spectral resolution of 1.5 A. As a rule, we observe two 
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spectral ranges: green, 4700-5400 A, containing strong absorption lines H/3, Mgib, FeIA5270 
and 5335, and red, 6200-6900 A, containing strong emission lines. Ha and [NII]A6583. The 
green spectral range is used to derive radial profiles of the absorption-line indices Mgb, 
Fe5270, Fe5335, and H/3 calculated in the well-known Lick system ([Worthey et al. 1994| ); to 



keep nearly constant signal-to-noise ratio along the radius, we have integrated the spectra 
in concentric circular rings with a width equal to one element (1'.'3) which are centered onto 
a galactic nucleus (for highly flattened systems, such as NGC 4216, this procedure results in 
smoothing details of the radial profiles). Also, we use the green spectral range to construct 
velocity fields of the stellar components; for this purpose we cross-correlate every galactic 
elementary spectrum with a spectrum of a bright star, typically K0-K3 III, obtained with 
the same instrumental setup. The emission lines in the red spectral range are measured 
to construct velocity fields of the ionized gas. We estimate the typical accuracy of the 
absorption-line indices to be 0.1 A and an intrinsic velocity accuracy ~ 20 km/s. 

To study in detail the morphology of the central parts of NGC 4216 and 4501, we have 
used a collection of digital images of 113 galaxies presented by Frei et al. (|1996|) . NGC 4216 



and 4501, among the others, were observed at the 1.5m telescope of Palomar Observatory 
on May 4, 1991. Three images, through the g, r, and i Thuan-Gunn filters, were obtained 
with the CCD TI of 800 x 800. The exposure times were taken equal to 60 sec, but the 
nuclei of both galaxies appeared to be saturated on the r- and z-images. The scale was 
1'.'19 per pixel, quite appropriate under moderate seeing conditions which accompanied the 
observations: Frei et al. ( |T996|) report seeing FWHM equal to 2'.'2-2'.'6 for NGC 4501 and 
2'.'6-3'.'4 for NGC 4216. Besides these moderate-resolution images, we have used an image 
of NGC 4501, obtained with the HST WFPC2 in the frame of the programme of Stiavelli 
"Core properties of the bulges of spiral galaxies", on May 3, 1997. The galaxy was exposed 
through the broad-band filter F606W, the total exposure time was 600 sec (400+200). 
The scale on the frame PCI which contains the nuclear part of the galaxy is 0'.'045 per 
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pixel, the spatial resolution estimated by a FWHM of foreground stars is O'/l, and under 
such sampling and resolution the central 35" x 35" of the galaxy has been exposed. The 
morphology of the images was analysed with the programme of one of us (V.V.V.) FITELL. 

3. Morphology of the Central Regions of NGC 4216 and NGC 4501 

3.1. Ionized gas 

Fig. 1 presents MPFS maps of the central parts of NGC 4216 and 4501 in the 
continuum and in the emission line [NII]A6583. It is clearly seen that in the center of 
NGC 4216 the emission-line morphology matches the continuum morphology, and the 
ionized gas distribution is similar to the distribution of stars. We agree with Keel (|1983|) 
who characterized the emission distribution in the center of NGC 4216 as a gaseous disk. 
Emission morphology in the center of NGC 4501 is less obvious, but we state that it is not 
point-like as Keel ( |1983| ) reported. We see an extended circumnuclear gaseous structure, 
but indeed the emission isophote shape in the center of NGC 4501 is more roundish than 
that of the continuum isophotes. Moreover, it lacks an axial symmetry: the northern part of 
this structure is fainter than the southern part. In other words, the observed morphology of 
the emission-line distribution in the center of NGC 4501 does not allow us to determine the 
spatial structure of the gaseous subsystem; it differs from the stellar distribution, although 
the kinematics of the ionized gas and stars in the center of NGC 4501 are quite similar (see 
the next Section). 

3.2. Stars 

Figs. 2-3 present some results of the isophote analysis using the central parts of the 
broad-band images of both galaxies. Here, together with our results, we plot K-band 
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isophote characteristics from the work of Rauscher (|1995|) - it is important for such dusty 
regions as the centers of NGC 4216 and 4501. In NGC 4216 there is no isophote major axis 
twist: a shght increase of PA toward the center according to the data of Frei et ah (|1996|) 



is not confirmed by the resuhs of Rauscher (|1995|) in the radius range of 1.5"-3". But 
other important features are evident: a local maximum of ellipticity at R ^ 11" and highly 
positive (particularly in the bluer filters) values of the fourth cosine Fourier coefficient 
at R <^ 11". The bulge subtracted photometric residual maps clearly show a highly inclined 
nuclear stellar disk with a radius of ~10"-11" which lies exactly in the plane of the main 



galactic disk. Rubin et al. ( |1997|) mentioned a weak bar in NGC 4216. But if it exists, 
it must lie outside the radius range under consideration; in the central part of NGC 4216 
no bar signatures are detected: edge-on bars look usually boxy, and in NGC 4216 within 
R ~ 30" there are no areas of negative a^. The edge of the nuclear stellar disk at i? ~ 10" 
can also be seen on the iJ-band surface brightness profile given for NGC 4216 in the work 
of Jungwiert et al. ( [L997|) . Jungwiert et al. (|1997|) looked for bars, either single or double, 
in disk galaxies; for this purpose they have made near-IR 2D photometry of 72 nearby 
galaxies and have analysed morphological characteristics of isophotes after deprojection of 
the images onto the galactic planes. NGC 4216 was not deprojected because of its large 
inclination, we could consider only the observed profiles. However NGC 4501 was done, and 
Jungwiert et al. ( |1997| ) concluded that in this galaxy there are no "bar or other triaxiality" , 
and slight oscillations of morphological characteristics near R ^ 18" are caused by spiral 



arms. Carollo et al. (|1998|) do not analyse their HST image of NGC 4501 because of 
multiple spiral arms penetrating into the very center of the galaxy. We have attempted to 
do this; together with the image characteristics derived from the Frei's et al. (|1996D data 
and the results of Rauscher ( |1995| ), the morphological characteristics of the isophotes of 
the HST WFPC2 image allow to make more detailed conclusions about the structure of 
the center of NGC 4501. The apparent ellipticity reaches a maximum value corresponding 
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to a thin disk inclined by 60°, namely, e=0.5, already at i? ~ 7", in the bulge-dominated 
region. Besides, the fourth cosine Fourier coefficient is always positive within R «i 9" 
according to the HST data and within R ^ 13" according to the data of Frei et al. (|1996|) 
and becomes negative outside this zone. Oscillations of PA plus-minus 6° from the line 
of nodes, seen even in the K-band, are more than probable within the zone of strong 04 
variations. So we detect morphological signatures of a circumnuclear stellar disk inside a 
radius of 10" and a slightly triaxial bulge dominating outside. As for nuclear "dust" spirals, 
they are clearly seen in the residual HST image of NGC 4501 obtained by subtracting a 
model image consisting of pure elliptical isophotes (Fig. 4). However, together with the 
"dark" spirals, "luminous" spirals are seen too. Perhaps, the SW, more luminous, arm 
causes the asymmetry visible in the emission-line image of the center of NGC 4501 (Fig. 1). 

4. Kinematics of the Central Regions of NGC 4216 and NGC 4501 

4.1. Stellar components 

Rotation of the stellar components in NGC 4216 and NGC 4501 has not yet been 
studied. We have obtained two-dimensional velocity fields of stars for the central regions of 
10" X 21" of both galaxies by using the green spectral range data from MPFS. They are 
shown in Fig. 5. 

First of all, two-dimensional velocity fields allow us to check if the rotation of stars is 
axisymmetric and if the plane of rotation coincides with the main plane of the galactic disk. 
When we have a disk which density (brightness) distribution and rotation are axisymmetric, 
the azimuthal dependence of the central line-of-sight velocity gradients must follow a cosine 
law within the solid-body rotation radius: 
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dvr/dr = uo sin i cos {PA — PAq), 

where uj is the deprojected angular rotation velocity, i is the inclination of the rotation 
plane, and PAq is the orientation of the line of nodes. Since the axisymmetric disk is 
round, the line of nodes is aligned with the isophotes. So if the maximum of the cosine-like 
azimuthal dependence of the central line-of-sight velocity gradients stays at the photometric 
major axis orientation, it strongly hints for an axisymmetric rotation. 

Fig. 6 presents the azimuthal dependence of the central line-of-sight velocity gradients 
for NGC 4216; it is cosine-like. Table 2 contains parameters of the best-fitting cosinusoid 
calculated by the least-square method. One can see that the orientation of the kinematical 
major axis is 201°. The orientation of the main major axis of the galaxy was known to be 
199° (RC3). We have obtained from the morphological analysis in the previous Section 
that the outermost isophotes have PAq = 201° and that the isophote orientation in the 
center is practically the same. So we conclude that the stars in the center of NGC 4216 
exhibit axisymmetric rotation. Taking into account that the galaxy is seen at an inclination 
of i = 78° ( [Guthrie 1992| ), we estimate the mean central angular rotation velocity as 28 
km/s/arcsec (or 350 km/s/kpc). The best circular-rotation model applied to the complete 
two-dimensional velocity field has given consistent values: PAq = 200° and i = 74°. 

Fig. 7 presents the azimuthal dependencies of the central line-of-sight velocity gradients 
for NGC 4501 according to the observational data from two runs, in May 1996 and 
May 1997. They are also cosine-like. Table 3 contains the parameters of the best-fitting 
cosinusoids. First of all, one can see that beyond R ~ 4" the cosine amplitude drops; that 
means that R ~ 4" is the border of the solid-body rotation region. Within R ^ 4" the 
mean value of PAq is 134° ± 3° which is not too different from the line-of-nodes orientation 
of 140° (RC3; our isophote analysis gave us 141°). Therefore we classify the rotation of 
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Table 1: Spectral observations of NGC 4216 and 4501 



Date 


Galaxy 


Configuration 


Exposure 


Field 


Sp. range 


PA 


Seeing 


12.05.96 


NGC 4501 


MPFS+CCD 520 x 580 


45 


min 


10" X 16" 


4800- 


-5400 A 


1° 


2" 


12.05.96 


NGC 4501 


MPFS+CCD 520 x 580 


40 


min 


10" X 16" 


6200- 


-6900 A 


1° 


2" 


15.07.96 


NGC 4501 


LS+CCD 520 X 580 


30 


min 


2" X 2.5' 


6200- 


-7000 A 


137° 


1.3" 


16.07.96 


NGC 4501 


LS+CCD 520 X 580 


33 


min 


2" X 2.5' 


6200- 


-7000 A 


46° 


1.3" 


16.07.96 


NGC 4501 


LS+CCD 520 X 580 


20 


min 


2" X 2.5' 


6200- 


-7000 A 


92° 


1.3" 


5.05.97 


NGC 4501 


MPFS+CCD 1040 X 1160 


60 


min 


10" X 21" 


4800- 


-5400 A 


1° 


2.5" 


2.05.97 


NGC 4216 


MPFS+CCD 1040 X 1160 


60 


min 


10" X 21" 


4800- 


-5400 A 


48° 


2.5" 


2.05.97 


NGC 4216 


MPFS+CCD 1040 X 1160 


40 


min 


10" X 21" 


6200- 


-6900 A 


48° 


2.5" 


21.01.98 


NGC 4501 


MPFS+CCD 520 x 580 


40 


min 


10" X 16" 


6200- 


-6900 A 


2° 


2.4" 



Table 2: Stellar rotation parameters for the center of NGC 4216 



Radius range, 


Amplitude, 


PAo, 


// 


cu sini, km/s/arcsec 


degrees 


1.5-2.4 


27.4 


201 
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Table 3: Stellar rotation parameters for the center of NGC 4501 



Radius range, 


Amplitude, 




// 


Lo sini, km/s/arcsec 


degrees 


May 96 






1.5-2.3 


29.6 


138 


2.7-4.2 


24.7 


128.5 


All 1.5-4.2 


26.2 


133.6 


May 97 






< 2 


28.3 


126 


2.6-3.7 


27.0 


138 


3.9-5.5 


19.5 


139 


All 1-3.7 


27.3 


134.1 
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the stellar component in the center of NGC 4501 as axisymmetric too. It implies that the 
angular velocity is 31 ± 1 km/s/arcsec (or 383 km/s/kpc) assuming i ^ 60°. 

Assuming that the rotation of stars in the centers of NGC 4216 and NGC 4501 is 
axisymmetric, we have calculated azimuthally averaged rotation curves in the frame of the 
circular rotation model. They are presented in Fig. 8. The parameters i and PAq were 
determined by fitting the model to the whole velocity fields. So for NGC 4501, the value of 
PAq is slightly different from that found by the cosinusoid analysis within the solid-body 
rotation area. However, this value of PAq is even closer to the main line of nodes. One can 
see that in Fig. 8 the region of solid-body rotation in NGC 4216 extends to R ^ 9", just the 
same as the radius of the circumnuclear gaseous disk (Rubin et al. |1997|); so the central part 



of NGC 4216 is kinematically decoupled from the more outer part of the galaxy. The radius 
of the solid-body rotation region in NGC 4501 is -R ~ 4" (Fig. 8). The local maximum 
rotation velocity of ~ 130 km/s at i? ~ 4" is followed by a slight velocity decrease. So we 
conclude that inside the radius of the kinematically distinct core in NGC 4501, 4" or 320 
pc, the mass is approximately 10^ Mq - a quite typical value for the kinematically distinct 
cores in spiral galaxies ([Afanasiev et al. 1989 ). 



4.2. Ionized gas 



The rotation of the ionized gas in NGC 4216 and NGC 4501 was studied more than 
once. We can mention the works of Rubin et al. ( [L989| ) and Woods et al. (|1990| ). These 
studies are fully based on long-slit spectra taken along the major axes; besides, they were 
intended to obtain large-scale, extended rotation velocity distributions, so their spatial 
resolution was rather low: Rubin et al. (|1989|) used a spatial binning of 1'.'5 per pixel, and 
Woods et al. (|1990|) observed under the seeing conditions of FWHM ^ 3" 5". Such an 
approach has allowed Rubin with co-authors to find a kinematically distinct circumnuclear 
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gaseous disk with a radius of 10" in NGC 4216 (Rubin et al. 1997 ), but a much more 



compact kinematically distinct core in NGC 4501 was not noticed. 

For the first time we present two-dimensional velocity fields of the ionized gas in the 
centers of NGC 4216 and 4501 (Fig. 9). If the gas velocity field for NGC 4501 looks quite 
regular and reminds us of the stellar component velocity field (Fig. 5b), the gas velocity field 
in the center of NGC 4216 is somewhat asymmetric. Taking into account that the galaxy 
is highly inclined and that we see the circumnuclear gaseous disk almost edge-on, we now 
wish to further study the gas rotation using the major-axis velocity measurements. For this 
purpose we have simulated " long-slit cross-section" by superposing a slit-like mask 1'.'3 x 20" 
on the two-dimensional velocity field. The results - one-dimensional velocity profile along 
the major axis - are presented in Fig. 10 together with the previous observations of Woods 
et al. ( |1990|) and of Rubin et al. ( |1997] ) (the systemic velocity of 130 km/s is subtracted 



from the measured line-of-sight velocities). Except for two points at R = 6" — 8" to the 
south-west from the center, the agreement is good. But in general, this line-of-sight velocity 
profile does not look regular. For the NE part the agreement is better, but Rubin et al. 



( [L997|) noted emission-line splitting at some radii. To the south-west from the nucleus the 
splitting was not noticed, but all three measurements diverge at i? > 2", implying that 
there may be some spatial sharp variations of the line-of-sight velocity pattern. Kuijken 
and Merrifield ( |1995| ) have simulated kinematical signatures of bars from edge-on to end-on; 
they have concluded that LOSVDs both of gas and stars (under the assumption of smooth 
spatial gas and star density distributions) must demonstrate gaps at some velocities giving 
a characteristic "figure-of-eight" picture. We have checked emission-line profiles (for the 
strongest emission line [NII]A6583) along the major axis of NGC 4216. With our spectral 
resolution, 1.5 A or 75 km/s, we clearly see LOSVD asymmetry with a long tail toward 
systemic velocity which is characteristic of the axisymmetric plane rotation seen almost 
edge-on; but we do not see any gaps in the profiles. So we conclude that there are no 
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kinematical signatures of a bar within R = 10" in NGC 4216 - this conclusion is consistent 
with the results of morphological analysis presented in the previous Section. The irregular 
shape of the major-axis velocity profile may be a consequence of a patchy dust distribution; 
some north-south asymmetry is also seen in the emission-line brightness map (Fig. la). An 
additional argument in favor of axisymmetric gas rotation in the center of NGC 4216 is the 
full agreement between the angular gas rotation velocity - 28 km/s/arcsec - reported by 
Rubin et al. ( |1997| ) (which is in fact line-of-sight velocity gradient along the major axis) and 
the stellar angular rotation velocity measured in the previous subsection; in the presence of 
a bar the stars (collisionless dynamic system) and the gas (dissipational dynamic system) 
should have different distributions and kinematics. 

Fig. 11 presents the azimuthal dependencies of central line-of-sight velocity gradients 
for the ionized gas of NGC 4501 in two different radius ranges and from measurements of 
various emission lines. The parameters of fitted cosinusoids are given in Table 4. First of 
all, one can see that the amplitude of the cosinusoids drops noticeably beyond R ^ 3": it 
is the border of the solid-body rotating gaseous disk. The data of May 1996 and January 
1998 are in agreement, and to determine parameters of the nuclear gas rotation, we involve 
measurements of [NII]A6583 in 1996 and of [NII]A6583, Ha, and [SII]A6717 in 1998. The 
orientation of the dynamical major axis for the ring of i? = 2"-3" is Py4o=141° ± 4°; it 
suggests that the rotation is axisymmetric, because the orientation of the main line of 
nodes is PAq = 140° (RC3) or 141° (our analysis). The averaged angular rotation velocity 
is 31 ± 1 km/s/" (or 388 ±17 km/s/kpc) if we assume the inclination to be 60°. We see 
that in the center of NGC 4501 both stars and ionized gas rotate axisymmetrically, in the 
equatorial plane of the galaxy, and with the same angular speed. 

Fig. 12 presents three long-slit cross-sections of NGC 4501: these profiles supplement 
major-axis cross-sections published earlier by Rubin et al. ( |1989| ) and by Woods et al. 
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Table 4: Gaseous rotation parameters for the center of NGC 4501 



Radius range, 


Emission line 


Amplitude, 








uj sini, km/s/arcsec 


degrees 


May 96 








2.0-2.4 


[Nil] 6583 


28.9 


147 


2.7-4.2 


[Nil] 6583 


20.4 


126 


Jan 98 








2.0-2.8 


[Nil] 6583 


24.2 


134 


2.0-2.8 


Ha 


25.7 


146 


2.0-2.8 


[SII]6717 


29.5 


134 


3.3-4.7 


[Nil] 6583 


16.4 


132 


All 
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( |1990|) . One can note that along the minor axis there is no velocity gradient near the center 
of NGC 4501; this is one more evidence for axisymmetric rotation of the circumnuclear 
ionized gas in NGC 4501. Fig. 13 shows the rotation curve for the ionized gas in the center 
of NGC 4501, both the azimuthally averaged one constructed from the two-dimensional 
MPFS data in the frame of an axisymmetric-rotation model and the one derived from the 
linear long-slit cross-sections. Up to i? ~ 3" the MPFS and long-slit rotation curves agree 
perfectly despite very different spatial resolution (seeing conditions) of two data sets, 1'.'3 
for the long-slit data and 2'.'4 for the MPFS data. But beyond R ~ 4" two branches of the 
long-slit major-axis cross-sections diverge extremely. By modelling the one-dimensional 
major-axis cross-section from the two-dimensional MPFS velocity field, we have obtained 
the same asymmetry. The azimuthally averaged MPFS rotation curve lies between these 
two branches, just as one could expect. The gaseous disk thus seems strongly asymmetric: 
perhaps, to the south-east from the center at i? > 3" we see the circumnuclear disk, and to 
the north-west - something else (dust screen?) that hides the NW part of the circumnuclear 
disk. If so, we must take into account only the SE branch of the major-axis cross-section, 
then the local maximum of the gas rotation velocity shifts to i? ~ 7", and the gaseous 
circumnuclear disk in NGC 4501 appears to be more extended that we thought before. 



5. Chemically Distinct Nuclei in NGC 4216 and 4501 

Figs. 14 and 15 present the radial profiles of the absorption-line indices for NGC 4216 
and NGC 4501, respectively. These profiles give us the evidence for a chemical distinctness 
of the nuclei in both galaxies: one can see a drop of metal-line equivalent widths when 
passing from the nuclei to the bulges. In the bulges which we take safely at i? > 4" (let us 
remind that seeing quality was 2'.'5) there are no noticeable gradients of metal-line indices, so 
we average the values in the radius range of 4"-9". Now we can make quantitative estimates 
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of differences in stellar populations properties by using models for old stellar populations. 
In NGC 4216 the mean metal-line indices in the bulge are < Mgb >= 4.20 ± 0.05 A and 
< Fe >= 2.72 ± 0.06 A (< Fe (Fe5270+Fe5335)/2); so the differences "nucleus-minus- 
bulge" are AMgb = 0.54 ± 0.15 A and A < Fe >= 0.74 ± 0.16 A. If we assume equal ages 
of the nuclear and bulge stellar populations and if magnesium-to-iron ratio is solar, then 
by using the models of Worthey ( P^^^ ) we obtain A[Fe/H]=0.26 ± 0.07 from AMgb and 
A[Fe/H]=0.41 ± 0.09 from A < Fe >. The difference between these estimates suggests a 
possible difference between magnesium-to-iron ratios in the nucleus and in the bulge. Before 
we make analogous calculations for NGC 4501, let us note that in the nucleus of NGC 4501 
{R < 3") a quite measurable emission line [NI]A5199 is seen in the spectra. As this line 
contaminates the red continuum passband of the Mgb index, its presence results in a Mgb 
overestimate. The effect has been discussed and quantified by Goudfrooij & Emsellem 
( |1996|) ; we apply their recommended correction, 1.13EW([NI])=0.45 A, to our nuclear Mgb 
measurement. After this correction, the bulge characteristics being < Mgb >= 4.26 ± 0.08 
A and < Fe5270 >= 3.08 ± 0.13 (according to the measurements of May 1997), we obtain 
the index differences AMgb = 0.56 ± 0.18 A and AFe5270 = 0.62 ± 0.23 A and, under the 
same assumptions as for NGC 4216, the metallicity difference 0.27 ± 0.09 from AMgb and 
0.38 ± 0.14 from AFe5270. Although the emission spectra of the nuclei of NGC 4216 and 
NGC 4501 are quite different, the properties of the stellar populations both in the nuclei 
and in the bulges appear to be very similar. 

But two main assumptions, that of equal ages and that of solar Mg-to-Fe ratio both 
in the nuclei and in the bulges, which allow to use the models of Worthey (|1994| ), are not 



necessarily valid. Let us try to clarify the situation in NGC 4216 and NGC 4501. For 
this purpose we can use "index-index" diagrams; but to refer to these diagrams, we need 
absolute values of the indices in the Lick system. Our spectral resolution in May 1997 
was 1.5 A - much better than that of the Lick system (8 A). So we may expect that our 
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index measurements in May 1997 are overestimated. Indeed, the differences between the 
measurements of NGC 4501 in May 1997 and in May 1996 have the right sense (Fig. 15). 
To check this effect, we observed bright stars (4 stars - in May 1997, resolution 1.5A, 9 
stars - later, resolution 4A) from the Lick library ( [Worthey et al. 1994|) . We have found 



that with the higher spectral resolution we overestimate the values of Mgb by 0.50 A and 
the values of Fe5270 by 0.44 A; with our standard resolution, which was used in May 
1996, the measured index values for the stars coincide with tabular ones within the errors 
of individual observations in the Lick project. Interestingly, just the same differences, 
0.4-0.5 A, can be seen between the May 1996 and the Mayl997 measurements in Fig. 15, 
so we conclude that the corrections cited above must only be applied to the measurements 
of May 1997 to transform them into the Lick system. 

Fig. 16 presents the diagrams (Fe5270, Mgb) for the central parts of NGC 4216 
and NGC 4501. One can see that these giant spiral galaxies differ strongly from a 
bulk of giant ellipticals on this diagram: if the latters lie to the right from the model 
sequences [Mg/Fe]=0 ([Worthey et al. 1992| ), the measurements for the nuclei and bulges 



of NGC 4216 and NGC 4501 lie close to them. The nucleus of NGC 4216 has obviously 
solar magnesium-to-iron ratio, the bulge may be slightly magnesium-overabundant. In 
NGC 4501, after applying the correction for the [NI] emission to the circumnuclear Mgb 
indices, both the nucleus and the bulge seem to lie close to Worthey's ( |1994| ) model 
sequences calculated for [Mg/Fe]=0. 

As the assumption of solar magnesium-to-iron ratio appears to be valid for the central 
parts of NGC 4216 and 4501, we can try to separate metallicity and age effects on the 
diagrams (H/?, Mgb), (H/?, [MgFe]) ([MgFe]= (Mgb < Fe >y/^), or (H/5, < Fe >). 
Unfortunately the absorption line H/? in the center of NGC 4501 is strongly contaminated 
by emission. But in the center of NGC 4216 the emission H(3 is practically absent - its 
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equivalent width is less than 0.1 A, - and we can estimate an age of the stellar population 
in the nucleus and its environment. Fig. 17 presents the comparison of our index data 
for the center of NGC 4216 with the models of Worthey (|1994|) on the diagrams (H/3, 



Mgb) and (H/5, [MgFe]) and with the models of Tantalo et al. (|1998D for [Mg/Fe]=0 on 
the diagram (H/5, < Fe >). Different metal absorption-line indices and different models 
give a stable result: the stars in the nucleus of NGC 4216 have a mean age of 8-12 billion 
years and a mean metallicity of 0.0 -r- +0.1, the stars in the bulge are older than 15-17 
billion years and more metal-poor than the Sun. So the chemically distinct nucleus in 
NGC 4216 is also younger than the nearest bulge. As for quantitative self-consistent 
estimates of the parameter differences, we can use recommendations of Tantalo et al. 
( p,998|) : they give linear equations which connect index differences AMgg, A < Fe >, 
and AH/5 with the stellar population parameter differences AlogZ, A[Mg/Fe], and logt. 
Assuming AMgb=16.2AMg2, we obtain for NGC 4216: A[Mg/Fe]=-0.12, AZ=+0.48, and 
Alogt=-0.39. If we take the direct Mg2 determinations from our spectra for the nucleus 
and for the bulge of NGC 4216 according to the prescriptions of Worthey et al. (1994), we 
obtain A[Mg/Fe]=-0.04, AZ=+0.49, and Alogt=-0.30, quite in agreement with the former 
results. So one can see that magnesium overabundance in the bulge is very modest, not 
more than 0.1 dex, the nucleus is younger than the bulge by a factor of 2-2.5, and, after 
taking into account the age difference, the nucleus is more metal-rich than the bulge by a 
factor of 3 instead of the factor of 2 obtained under the assumption of equal ages. 



6. Conclusions 

By using bidimensional spectral data obtained at the 6m telescope for the Virgo 
spirals NGC 4216 and NGC 4501, we have found chemically distinct metal- rich nuclei 
in these galaxies. Under the assumption of equal ages for the nuclear and bulge stellar 
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populations, the metallicity difference between the nuclei and their environments in the 
galaxies is estimated as a factor of 2. But we have also found an age difference between 
the nucleus and the bulge in NGC 4216: age-metallicity disentangling results in an 
age estimate for the nucleus of 8-12 billion years, the bulge being older by a factor of 
1.5-2; and the self-consistent metallicity difference estimate is then a factor of 3. The 
solar magnesium-to-iron ratios in the galactic nuclei suggest for a long duration of the 
secondary nuclear star formation bursts which have produced the chemically distinct stellar 
subsystems. 

The detailed morphological and kinematical analyses made for the stellar and gaseous 
structures in the centers of NGC 4216 and 4501 have revealed the presence of the 
circumnuclear stellar-gaseous disks which exhibit fast axisymmetric rotation and lie in their 
respective equatorial planes of the main galactic disks. The radius of the circumnuclear 
disk in NGC 4216 is 10" (800 pc), and this value agrees with the estimate of Rubin et al. 
( |1997|) made for the purely gaseous inner disk. The inner disk radius in NGC 4501 can only 
be roughly estimated in the range 4"-7" due to the complex dust distribution in the central 
part of the galaxy. Unlike the elliptical galaxies with chemically distinct cores where area 
of magnesium-index variations coincides with the inner-disk extensions ( [Surma fc Bender 



19951 ), the chemically distinct nuclei in NGC 4216 and 4501 seem to be unresolved although 
the inner disks have quite measurable extensions. This work presents the second evidence of 
this sort; the first was the case of M 31 (ISiPchenko et al. 1998| ), but it looked more striking 



- a circumnuclear 100 pc stellar disk in this galaxy extends up to 0.'5 from the center, and 
the chemically distinct nucleus is much more compact, less than 4". One may note that 
the inner disks contribute less than 20% of the total light in the bulge-dominated areas so 
they cannot affect significantly the integrated spectra of the stellar populations. But in 
the elliptical galaxies the situation is similar, meanwhile the disks are detectable by their 
magnesium-enhanced contributions. In any case, the unresolved chemically distinct nuclei 
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in spiral galaxies represent something different from the inner disks, though there may exist 
some genetic relation between these two substructures. 
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Fig. 1. — Surface brightness distributions of the continuum at A6500 (grey-scaled) and the 
emission [NII]A6583 (isophotes) for the central parts of NGC 4216 (a) and NGC 4501 (6) 

Fig. 2. — Results of the isophote analysis for the broad-band images of the central part of 
NGC 4216: a — major axis position angle radial variations, h — ellipticity radial variations, 
c — boxiness radial variations. Asymptotic parameter values R — 120" — 200" are: 
PAq = 21°, e = 1 - 6/o = 0.83 

Fig. 3. — Results of the isophote analysis for the broad-band images of the central part of 
NGC 4501: a — major axis position angle radial variations, h — ellipticity radial variations, 
c — boxiness radial variations. Asymptotic parameter values at i? = 120" — 160" are: 
PAq = 141°, e = l-h/a^ 0.54 

Fig. 4. — The residual map (after pure elliptical isophote subtraction) of the HST WFPC2 
image of NGC 4501. The hght ellipse marks the region where model subtraction has been 
made; the image is direct, PA{top) — 15?3, the total image sizes are 36" x 36". Luminous 

(grey) and dark (white) spiral arms are clearly seen in the radius range of 1" — 10" 

Fig. 5. — Two-dimensional line-of-sight velocity maps for the stellar component in the centers 
of NGC 4216 (a) and NGC 4501 (6). The images are direct, the north is pointed by an arrow, 
the photometric center is marked by a cross 

Fig. 6. — Azimuthal dependence of the central stellar line-of-sight velocity gradient for 
NGC 4216 

Fig. 7. — Azimuthal dependence of the central stellar line-of-sight velocity gradient for 
NGC 4501, in five radial bins 

Fig. 8. — Azimuthally-averaged rotation curves for the stars in the centers of NGC 4216 
and 4501 calculated over the observed 2D MPFS velocity fields in the frame of a two- 
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dimensional axisymmetric model; the best-fit parameters for the orientation of the galaxies 
are PAq = 200°, i = 74° for NGC 4216 (a) and PA^ = 140°, i = 61° for NGC 4501 (b). For 
comparison the major-axis velocity profiles simulated by masking our 2D velocity fields are 
also plotted as sohd hnes 

Fig. 9. — Two-dimensional line-of-sight velocity maps for the ionized gas in the centers of 
NGC 4216 (a) and NGC 4501 (6). The images are direct, the north is pointed by an arrow, 
the photometric center is marked by a cross 

Fig. 10. — The major-axis hne-of-sight velocity profile for the ionized gas ([NII]A6583 
emission line) in the center of NGC 4216 simulated by masking our 2D velocity map. For 
comparison the long-sht observations of Woods et al. (1990) and of Rubin et al. (1997) are 
also plotted 

Fig. 11. — Azimuthal dependencies of the central line-of-sight velocity gradient for the 
ionized gas in NGC 4501 in several radial bins; sohd lines show the best-fit cosinusoids with 
parameters listed in Table 4 

Fig. 12. — Long-slit profiles of ionized-gas hne-of-sight velocities in PA — 137° (a), PA — 46° 
(6), and PA = 92° (c) for NGC 4501 

Fig. 13. — The ionized-gas rotation curve for NGC 4501. The points with error bars are 
calculated from the 2D MPFS velocity fields in the frame of a two-dimensional axisymmetric 
model; the best fit is reached for the orientation mentioned in the legend. The long-sht 
major-axis measurements divided by sin i are also plotted; a strong discrepancy between the 
SE and NW parts of the major-axis cross-section is observed ai R> 3" 

Fig. 14. — Radial profiles of the absorption-line indices in NGC 4216 

Fig. 15. — Radial profiles of the absorption-hne indices in NGC 4501 
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Fig. 16.— The diagram (Fe5270, Mgb) for the central parts of NGC 4216 (a) and NGC 4501 
(b). The observational points are taken along the radius with a step of 1'.'3 and connected 
by dashed lines. The ages of the model stellar populations are given in the legend in billion 
years 

Fig. 17.— The diagrams (H/?, Mgb) (a), (H/?, [MgFe]) (6), and (H/3, < Fe >) (c) for the 
central part of NGC 4216. The observational points are taken along the radius with a step 
of 1''3 and connected by dashed hues. The ages of the model stellar populations from the 
works of Worthey (1994) {a,b) and Tantalo et al. (1998) (c) are given in the legend in billion 
years 
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